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Abstract 

We apply the recently defined Lambert W function to some problems of classical statistical me- 
chanics, i.e. the Tonks gas and a fluid of classical particles interacting via repulsive pair potentials. 
The latter case is considered both from the point of view of the standard theory of liquids and in 
the framework of a field theoretical description. Some new mathematical properties of the Lambert 
W function are obtained by passing. 



PACS numbers: 05.20.-y, 64.10.+h 



I. INTRODUCTION 



The Lambert W function is defined as the multivalued inverse of the function w i— > we w . 
Its mathematical properties have been explored only quite recently after its implementa- 
tion in the mathematical library of the computer algebra program Maple. The history, 
applications and properties of Lambert W are reviewed in Corless et a/.— Mathematical de- 
velopments and applications to physics (mainly in quantum mechanics and electrostatics) 
can be found in refs. Q^l^U. Since W is a very simple function, all its applications to physics 
or other fields are not exhausted, and, after a short enquiry, I have discovered that many of 
my colleagues, aware or unaware of its name, have met the function in their own works. 5 

Here we reformulate some old and new problems of the statistical mechanics of classical 
liquids in terms of Lambert W. 

Our paper is organized as follows. In sec. II we commence by refreshing the reader 
ideas with some known mathematical properties of W(z). We focus on the principal branch 
Wq(z) of the Lambert function and the related function Uq(z) = Wq(z) + Wq(z)/2, both 
being of prime importance for their applications to the problems of statistical mechanics 
considered in this article. Some new mathematical properties of W (z) and U (z) used in 
the paper are presented in the appendix (the Legendre transforms of the real functions 
x G E ,x i— > Wo(exp(x)) and x i— ► Uo(exp(x)) are computed and the dispersion relations for 
W Q (z) and U (z) established). 

In sec. Ill we reformulate the properties of the Tonks gas,— i.e. a one- dimensional (ID) 
classical fluid of hard rods, in terms of the principal branch Wq{z) of the Lambert function. 
Our presentation is compared with the seminal work of Hauge and Hemmer on this system.— 
The introduction of Lambert Wq(z) allows to recover all the properties of the model in a 
unified treatment. 

In sec. IV we focus on the second system considered By Hauge and Hemmer in ref. [D, i.e. 
a classical 3D fluid made of particles with intermolecular pair repulsions of the form 

<f(r) = 1 3 0M (1.1) 

where is a positive integrable function. The extension to arbitrary dimension is trivial. In 
the limit 7 — > which will be considered, we deal with infinitely weak and infinitely long 
range repulsive interactions. Once again, the use of the functions Wq(z) and Uq(z) is of a 
great help for describing in a unified manner all the known properties of this model. 



Sec. V contains new material. We consider the field theoretical representation of a fluid 
of particles interacting via repulsive potentials in the framework of the formalism developed 
in ref. It is shown that at the mean field (MF) level of the theory the expressions for the 
pressure and the density are identical to those derived in sec. IV for the fluid with weak long 
range repulsive interactions. In other words, the MF approximation is exact in the case of 
infinitely long range interactions, a satisfactory, although expected result. We conclude in 
Sec. VI. 



II. A DIGEST ON LAMBERT W 



In this section we summarize the main mathematical properties of the Lambert W(z) 
function. W(z) is defined as the root of 

W{z)e w{z) = z , (2.1) 



where z is a complex number. W(z) is a multivalued function which has been studied 
recently in Refs. The different branches of Lambert W (i.e. the different possible 

solutions of eq. (J2.1)0 are labeled by an integer k = 0, ±1, ±2, etc. When z is a real number 
eq. 1)2.1 j) can have either two real solutions for > z > — exp(— 1), in which case they are 
W / o(^) an d W-i(z), or it can have only one real solution for z > 0, this being Wq(z) while 
W-i(z) is now complex, or no real solutions for — oo < z < — exp(— 1). Wq{z) and W-i(z) 
are the only branches of W that take on real values. 

A brief survey of the properties of the principal branch Wq(z) will be sufficient for our 
purpose. Firstly, Wq(z) is analytic at z = 0. This follows from the Lagrange inversion 
theorem^. Its power series about the point z = reads as 

W (z) = ^^—^. (2.2) 

' 71! 

n=l 

The radius of convergence of the series (j2.2j) is equal to e~ l as it is easily shown using the 
Cauchy test and, within the circle of convergence, |W(z)| < 1. As a consequence of the 
relation (|2.1j) one has 

w »w = i(Mm • (2 - 3a) 



< 2 ' 3b > 



> 1) 

n=l y 



where the radius of convergence of the series (|2.3b|) is again equal to e -1 . Lagrange theo- 
rem gives more than eq. ([2.2)1 : it allows to obtain power series for Wq(z), 1/(1 + Wq(z)), 
exp(aWo(z)), etc for z within the circle of convergence of Wo(z)£ In particular, the power 
series of the function Uq(z) defined to be Wq(z) + Wq(z) 2 /2 is given by 



U (z) = { ) , (M < e" 1 ) • (2.4) 

71=1 

Note that if follows from eq. (|2.3a|) that 

zU' Q {z) = W (z) . (2.5) 

Wq(z) has a second-order branch point at z — — e _1 which it shares with both W\{z) and 
W-i(z) and its branch cut is conveniently chosen to be {z : — oo < z < — e -1 } with the 
convention that Wq(z) is defined on the upper lip of the branch cut. The behavior of Wq(z) 
about the branch point is given by the series^ 6 

W (z) = -l+p-^p 2 + ^p 3 + ... , (2.6) 

where p = a/ (2(ez + 1) (the series converges for \p\ < \/2). 

In physics, it is quite common (and sometimes enlighting) to write dispersion relations 
for functions having a branch cut such that Wo(z) and Uq(z)^ We show in the appendix 
that, in the present case, these relations can be recast under the form 

r e_1 z 

W(U) Q (z) = / g w{u) (s)\n(l-~)ds, (2.7a) 

9w(u)(s) = -1±Q\W{U){ 8 )] , (2.7b) 

where z is an arbitrary complex number which however is not on the branch cut. These 
equations convey the interpretation of Wo(z) (resp. Uq(z)) as the two-dimensional complex 
electric potential created by a distribution of charges gw( s ) (resp. gu( s )) located on the 
branch cut. Moreover, we show in the appendix that the distribution gw{ s ) is normalized to 
unity (cf eq. ([A. 14)) ) while the total charge of the charge distribution gu{s) is infinite. Other 
physical interpretations of the relations ()2.7j) will be discussed in forthcoming developments. 



III. THE TONKS GAS 



In this section we consider the Tonks g ID fluid of hard rods of length <r.— The 

equation of state (EOS) is known exactly: 

X = /?P = — . (3.1) 
1 — pa 

(P pressure, j3 = 1/kT, T temperature, p number density). Of course < pa < 1. For 
pa — > one recovers the EOS of the ideal gas, while for pa — > 1, which corresponds to the 
close packing of the hard rods, the pressure diverges. Less trivial is the expression of the 
pressure as a function of the activity z. Recall that z = exp(/3/i) where p is the chemical 
potential and we have assumed that the deBroglie thermal wavelength A = 1. Starting from 
the thermodynamic relation 

P = z j-Xi ( 3 - 2 ) 

it is not difficult to obtain the relation 7 

a X e ax = az . (3.3) 

A priori Eq. ()3.3|) is valid only for real positive values of z but it allows to define the pressure 
x{ z ) i n the complex plane of the activities by analytic continuation. The authors of ref. Q 
did not know the Lambert function as we do; clearly one has 

X (z) = -Wo(az) , (3.4) 
a 

since the other branches of Lambert W will not give a real pressure for real chemical 
potentials p. The Mayer expansions for the pressure and the density of the Tonks gas are 
therefore obtained by copying out eqs. ()2.2|) and (j2.3b|) : 



oc 



X (z) = Y,Kz n , (3.5a) 

n=l 

1 W (az) 
a 1 + W (az) 



1 W Q (az) ^ 
p(z) = - ^ WJnA =l^nb n z , (3.5b) 

n=l 

(-na) n - 1 

b n = r — ■ 3 - 5c 

n! 

The results of sec. II on Lambert W enable us to conclude that 

1. The radius of convergence of the Mayer series ()3.5a|) and (|3.5b|) is R = 1/ea. Within 
the circle of convergence \x( z )\ < & 1 - 



2. The pressure x( z ) is singular at the point z = —R of the circle of convergence, i.e. 
the branch point of Wo(az), in agreement with one of the conclusions of the second 
theorem of Groeneveld.— 

3. x{ z ) has a branch cut on the negative real axis: — oo < s < —R. It can be identi- 
fied with the distribution of zeros of the grand partition function in the thermody- 
namic limit . 12 i 13 The distribution of the Yang-Lee zeros is given by g(s) = gw{°~ s ) — 
— ^sW^as) /n (cf eqs. (I2.7a|) and (jA.13a|0 . g(s) is a positive and increasing function 
over the interval — oo < s < —R which behaves as g(s) ~ — l/(asln |s|) 2 ) for s — > — oo 
and as g(s) ~ ^ e/2/(irx 1 / 2 ) for x = —R — x, x — > 0+; therefore the integral of g(s) 
over the interval — oo < s < —R is convergent and equal to a (cf eq. (jA.14p ). 

4. For an arbitrary z (however not on the branch cut) the pressure xi z ) may be seen 
as the 2D complex electrostatic potential created by the charge distribution g(s). It 
follows from eqs. ()2.7a|) and (jA.13a|) that x{ z ) takes the form proposed by Yang and 
Leei^ 

r R z 

x(z) = / g(s) ln(l ) ds (Vz not on the cut — oo < s < —R ) (3.6) 

J -oo S 

A close examination of their paper reveals that the derivation of eq. ()3.6J) given by 
Hauge and Hemmer is restricted to a point z inside the circle of convergence. 

It is amusing to check the various conclusions of second theorem of Groeneveld.— Recall 
that for a d-dimensional classical fluid of particles interacting via positive pair potentials 
fi'f'ij) ^ such that (twice) the second virial coefficient 

f = J rf d f(exp(-Mr))-l), (3.7) 

converges, the radius of convergence R of the virial series (J3.5a|) and the Mayer coefficients 
b n satisfy to the following inequalities 



^<*<i (3.8a) 

n j n - x ~ n\ v ' 

In the present case / = 2&2 = — c, and one easily checks that R = 1/ea and b n as given by 
eq. (|3.5c|h do satisfy to the Groeneveld inequalities ()3.8|) . 



We end this section by determining the Helmoltz free energy per unit volume f3f{p) 
as the Legendre transform of /3P viewed as a function of the reduced chemical potential 
v = (5p = Inz (see eg ref. I14L It follows from the eq. ()A.5|) of the appendix that 

Pf(p) = sup{vp- X {e u )) 

i/<EK 

= p(lnpA- 1) -pln(l -pa) . (3.9) 

Note that we have restored A to make the argument of the In dimensionless. (3f(p) is a 
strictly convex function of the density defined for < per < 1. The first term in the RHS of 
the equation is the Helmoltz free energy per unit volume of the ideal gas while the second 
one is the excess free energy. The analytic continuation of f3f{p) to complex densities is 
obtained by defining f3f(p) according to eq. ()3.9j) where In is the principal branch of the 
natural logarithm for (3f(p) must take on real values for p £ M, < pa < 1. f3f{p) is then 
a multivalued function with two branch cuts —00 < p < and a^ 1 < p < +00. 



IV. WEAK AND LONG RANGE REPULSION 

The EOS of a gas of particles interacting via the pair potential (jl.lj) is known exactly 
and given by£ 

X = /3P = P+fp 2 , (4.1) 

where 

a = /?0(O) = J d 3 rf30(r) (4.2) 

is a positive, 7-independent constant. Here p can take on all non negative real values. It is 
not difficult to obtain the relation between the density and the activity which reads as^ 

z = pe ap . (4.3) 
Therefore, since for a real z the density must be real 

p = -W (az) , (4.4) 
a 

and, by integration of the thermodynamic relation ()3.2|) 

X {z) = -Uo(az) . (4.5) 
d 



The pressure and the density for complex activities are obtained by an analytical continu- 
ations of eqs. (|4.4|) and ()4.5|) . The virial series for x{ z ) an d Pi z ) follow from eqs. (|2.2|) and 
(El 



-a) n - l n n - 2 



X (z) = £ ^ ' , z» , (4.6a) 

n=l 

?w = E — „, ■ <« b > 

n=l 

The radius of convergence of these two series is R = (ae) _1 and it satisfies to the first 
Groeneveld inequality (|3.8a|) . One also checks that the coefficients of the power series (I4.6a|) 
do satisfy to the second Groeneveld inequality (|3.8b[) . x( z ) nas a branch cut on the negative 
real axis: — 00 < s < —R. For any complex activity z not on the cut, it can be written once 
again under the form proposed by Yang and Le o 12 ' 13 (cf eq. (jA.13b)0 

X (z)= [ R 9(s) ln(l - -) ds , (4.7) 



where 



g(s) = gu (as) = -—%W (as) . (4.8) 



It follows from the analysis given in the appendix that the function g(s) is not integrable on 
the cut (it behaves as —1/as for s — > —00). Recall that, strictly speaking, the theory of Yang 
and Lee is valid only for interactions with a hard core contribution. In this case the grand 
partition function 5 for a finite volume V is a polynomial in z. The distribution g(s) of its 
zeros after the passage to the thermodynamic limit allows to rewrite x{ z ) = InH/V under 
the form (|4.7j) . where g(s) is normalizable . 7 i 12 i 13 For soft interactions as those considered 
in this section, the reasoning breaks down, leaving us with a non normalizable distribution 
g(s) (one can pack an infinite number of particles in a finite volume and the density p is not 
bounded). 

Finally the Helmoltz free energy per unit volume (3f(p) is computed as the Legendre 
transform of x( e ^)- ft follows from the eq. (|A.6|) of the appendix that 

Pf(j>) = sup(up--U (aen) 

= p(lnpA 3 -l) + ^p 2 . (4.9) 

Pf(p) is a strictly convex function of p on the interval < p < +00. 



V. REPULSIVE INTERACTIONS: A FIELD THEORETICAL APPROACH 



We consider the statistical mechanics of a system made of N classical point particles 
interacting via pair potentials of the form 

v(r) = <p (r) + <p(r) , (5.1) 

where (fo(r) is some reference potential (eg a hard sphere repulsion for instance) and (f(r) 
admits a positive, well defined Fourier transform <p(k) > 0. <po(r) an d <f(r) are supposed to 
meet all the requirements which are necessary for the existence of a thermodynamic limit 
(TL)M 

We denote by Q the domain occupied by the molecules of the fluid. It will be convenient 
to assume that Q is a cube of side L with periodic boundary conditions (PBC). In a given 
configuration u = (N; r i . . . r N ) the microscopic density of particles reads as 

JV 

p{f ) = Y^8^{r-r i ) , (5.2) 



and its Fourier transform 



i=l 



N 



Pk = y^exp(-ifc.f i ) . (5.3) 



i=l 



The configurational potential energy of the system can be decomposed as 

PV(u) = PV (U) - Nu s + |<pMp) , (5.4) 

where Vo(uj) is the configurational energy of the reference system, v$ = /3<p(0) a self-energy 
contribution and 



(p\(p\p) = / d rid f 2 p(rl)^(fi,f 2 )p(rl) . (5.5) 
Jn 

We shall work in the grand canonical (GC) ensemble. We denote by [i the chemical po- 
tential and by ip(f) an external potential. The local chemical potential will be noted 
v{f) = /3(fi — ip(r)). Performing a Kac-Siegert-Stratonovich-Hubbard-Edwards (KSSHE) 
transfor m 16 i 17 i 18 i 19 i 20 one can show that the GC partition function can be recast under the 
form of a functional integral 8 

~\v]=H- 1 j Vtexp{-H[{}), (5.6) 



where X>£ is the functional measure and the action of the KSSHE field theory reads as 



m = h ~ l0gH ° [V+ iC] ' (5 ' 7) 

where V = v + vg and So is the GC partition function of the reference system. The field £ 

which enter eqs. (|5.fi|) and (|5.7|) is a real scalar field. Note that, in the general case where the 

sign of (p(k) is abitrary, the action TC involves two real scalar fields and £_ (or a complex 

field £).— Finally, the normalization constant M is given by 

AT= J vtexpi—itlip- 1 ®). (5.8) 

In ref. the potential (po(r) of the reference system was chosen to be a hard core potential 
of diameter a. Here we specialize to the case a — > (or equivalentlv 9? = 0), i.e. we take 
for the reference system the ideal gas. All the conclusions of ref. |s| remain valid provided 
that the TL of our system is well behaved. Note that the configuration energy (|5.4|) can be 
rewritten as 

nu) = ^^\ps\ 2 m - %<p(p) - M) 

Therefore if (p(k) > for all k and ip(0) > (which will be assumed henceforth) then 
V(oS) > —NB with B{= </?(0)/2) > 0, i.e. the system is H-stable in the sense of Ruelle 
and the TL exists.— Note that, conversely, if <p(k = 0) < 0, the system does not have 
a thermodynamic behavior and the introduction of a repulsive hard core is mandatory to 
ensure the existence of a TL. Such potentials (Le. hard core plus an attractive tail) can yield 
a liquid-vapor transition as explained in ref. |8j in the case considered here (soft repulsive 
tail), the possibility of such a transition has to be ruled out. 
With the choice (po = the KSSHE action reads now 

= 4 (^-'10 ~ [ d 3 fe v ^m . (5.10) 

We turn now our attention to the mean field (MF) level of the theory. The MF or saddle 
point approximation is defined by the equation 

Z MF (v) = exp(-W(0) , (5.11) 

where at ^ = ^ the action 7i is minimum. If there are several local minima then one retains 
the absolute minimum. The stationary condition 

5H 



(5.12) 



can be recast under the form of the implicit integral equation 

f (r) = ip f <p(\\r - f*\\) e 17 ^^ . (5.13) 



Jn 

Moreover the MF density of the fluid is given by the density of the reference fluid -here the 
ideal gas- at the local chemical potential V(f) + i£(r)£, i.e. 

p MF (f) = e^+^W . (5.14) 

In the case of a homogeneous system to which we stick from now, (i.e. ip(r) = 0) the 
solution of eq. ()5.13|) is clearly 

£ = iW k (\z), (5.15) 
where Wk is some branch of Lambert W, z the activity and 

A = pfi(0)e us , (5.16) 

from which it follows that 

Pmf = W k (\)//3<p(0) . (5.17) 

It remains to determine the branch of W. For real activities z, the MF density should be 
a positive real number; it follows then from eq. (J5.17j) that, necessarily, k = 0. The MF 
pressure is easily derived from eq. (|5.11|) with the result 

Xmf{z) = \nE MF /Q , 
1 



Note that the MF equation of state takes the familiar form Xmf(p) — P+/3^p{0)p 2 /2. Finally, 
the MF free energy is obtained by a Legendre transform of x( z )> which gives 

(3f MF = p(ln(pA 3 ) - 1) - pu s + ^(0)p 2 /2 , (5.19) 

in agreement with the general expression of ref. 0. 
Some comments are in order. 

1. The Taylor series of xmf(z) and Pmf(z) in terms of the activity z are once again 
given by eqs. (|2.2jl and ()2.4j) with a radius of convergence equal to R = (Ae) -1 . For an 
arbitrary complex activity z, the pressure Xmf( z ) is given by the dispersion relation 



2. We have shown elsewhere^ that the one-loop order approximation of the KSSHE 
field theory coincides with the random phase approximation (RPA) of the theory 
of liquids. 21 As well known, for long range potentials of the form the RPA cor- 
rections to the pressure vanishes as 7 — > O. 21 Note that, in this limit, z/g — * and 
£>(0) = 0(0); thence, in the limit 7 — > 0, one recovers for the pressure, the free energy 
etc the expressions derived in sec. IV. Stated otherwise, the MF-KSSHE theory is 
exact for infinitely weak and long range repulsive potentials. 

3. We have shown in ref.|8|that the MF pressure and free energy constitute exact bounds. 
More precisely we have 

\fzER X {z) < Xmf{z) = U (\z)/P^(0) , (5.20a) 
Vp > (3f(p) > (3f MF (p) = p(ln(pA 3 ) - 1) - pu s + /^(0)p 2 /2 . (5.20b) 



VI. CONCLUSION 



In this paper we have discussed some applications of the Lambert W function to the theory 
of liquids. In the case of ID hard rods or infinitely weak long range repulsive potentials, one 
obtains a close expression for the complex pressure x( z ) as a function of complex activities 
in terms of either Wq(z) or the related function Uq(z). The dispersion relations derived 
for Wq (z) and Uq(z) in the appendix give a rigorous justification to the heuristic formula 
proposed by Lee and Yang for x(^) . 12 ' 1 ? We have also shown that, in the framework of the 
KSSHE field theory of liquids, the MF pressure x( z ) °f a g as °f particles interacting via soft 
repulsive potentials can be expressed in terms of the function Uq(z). 
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APPENDIX: NEW MATHEMATICAL PROPERTIES OF LAMBERT W 



1. Legendre transforms 

Let us first consider the real function x \— > f(x) = W (e x ). By applying twice the formula 
(|2.3a|) one finds 

/'M - ^ , (A.la) 

f» - (rrf)F • (A - lb) 



It follows from eq. (jA.lb|) that f"(x) > for all x G R. Thence the 



unction / is strictly 



convex and one can define its Legendre transform Cf(p) (see eg ref. Il4). By definition 

£ / (p)=sup(xp-/(x)) > (A.2) 

or, more precisely 

C f (j>) = xp-f{x), (A.3) 

where x is (the unique) solution of 

Therefore Cf{p) is defined for < p < 1. It is also a strictly convex function, the Legendre 
transform of which is f(x) (i.e. the Legendre transform is involutive). In order to get 
the expression of £f(p) one notes that eq. (jA.4|) implies that fix) = p/(l — p) and that 
x = fix) + Infix) which yields 

C f {p) = pOnp-l)-phi{p-l). (A.5) 

Similarly it can easily be shown that the function x i— * h{x) = Uq{c x ) is strictly convex 
and that its Legendre transform Chip) * s given by 

C h (p)=pQnp-l) + £. (A.6) 

Note that Chip) * s defined on < p < oo). Over this interval, Chip) * s strictly convex. 



2. Dispersion relations 



We establish here the dispersion relations for the functions Wq(z) and Uq(z). Let us first 
consider Wq(z). It follows from the Cauchy theorem that 



where C is the contour shown in fig. (JIJ) and z is not on the cut (— oo, — e _1 ). 
Wo(0) = as follows from eq. (|2.2j) . We consider now eq. ()A.7|) in the limit e - 



(A.7) 

Indeed 
> and 



oo. 




FIG. 1: Solid line: contour of integration of eq. (|A.7|) . The small circle (dashed line) is the circle 
of convergence of Wq (z) about z = 0, the dot is the branch point and the thick solid line is the 
branch cut. The point z is everywhere except on the cut. 



The asymptotic formula for large (complex) z, i.e 



1,4 



W {z) ~ ln(z) - ln(ln(z)) 



(A.8) 



where In z is the principal branch of the natural logarithm, ensures that the contribution 
to eq. (jA.7|) from the large circle T tends to zero as its radius R tends to infinity (Jordan 
lemma). Note that we have substracted Wo(0) = from W / o(^) in eq. (|A.7|) precisely in 
order to obtain this property. Therefore we have 



Wo{z) = ^-. [ £ { W (s + ie) ' ' 1 



s + ie — z s + ie 



- W (s - ie) [ — ] } ds . (A.9) 

s — it — z s — ie 

As e — > 0, we can neglect the ±ze in the fractions that enter the RHS of eq. (|A.9J) (remember 
that z is not on the cut) and we get 

Wo(z) = -^ f [Wo(s + ie)-Wo(s-ie)][—--}ds. (A.10) 
2m y_ 00 s-z s 

With the convention that Wo( s ) is defined on the upper lip of the cut and by noting that 
Wo(s) = Wo(s) (for s not on the cut) we infer from eq. (jA.10|) that 

W (z) = -[ £ Q(W { 8 ))[—--]d8. (A.ll) 

7T S-Z S 

The final step is to perform an integration by parts which yields 



W (z) = [3(Wo(s))ln(l--)] 



Z^-e- 1 
oo 



S 

~ e ~' -H„ z 



— oo 



. »(W (s)) ln(l - -) ds . (A.12) 
7i ds s 



Since, in the one hand Wq{— e" 1 ) = —1 and, on the other hand, 3(Wo(— oo)) = vr as 
can be obtained for instance from the asymptotic behavior (jA.8|) of Wq(z), then the first 
contribution to the RHS of eq. (|A.12|) vanishes and we are left with 

r e_1 z 

W (z) = / g w {s) ln(l - -) ds , (A.13a) 

J —oo S 

9w(s) = -~*\Wo(8)] ■ (A.13b) 
Ti as 

Several comments are in place here. 

(a) The function gw( s ) is a positive increasing function on the cut (— oo, — e _1 ). For 
s = — e _1 — x (x — > +0), gw(s) ~ \/e/2/ (ttx 1 ^ 2 ) as can be seen from the series expansion of 
Wo(z) about the branch point (cf eq. (|2.fij) ). In the other hand gw( s ) ~ — l/(s In |s|) 2 ) for 
s — > — oo as can be infered from eq. (|A.8j) . Thence the integral of gw( s ) along the branch 
cut is convergent; more precisely 

/ g w (s) ds = [Q (W (-oo)) - 3 {Wo (-e" 1 ))] /tt 

J — oo 

= 1 . (A.14) 



(b) The dispersion relations for the function Uq(z) are similar to eqs. (jA.13aD for the 
reasoning presented above can be reproduced without major changes. The distribution 
gu{s) is again a positive function on the interval (—00, — e _1 ). It is given by 



gu(s) = -7T 1 d(QU (s))/ds 

= -95«i> , (A.15) 

TTS 

from which it follows that gu{~ el ) — because W^—e^ 1 ) = — 1. Moreover, for s — ► —00 
we have ~ — 1/s with the consequence that the integral of gu{s) along the branch cut 

is divergent to +00. 
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